The Extragalactic Infrared Background and its Cosmological Implications 

IAU Symposium, Vol. 204, 2001 

Martin Harwit and Michael G. Hauser, eds. 



Global Evolution Models and the Extragalactic 
Background Light 



S. Michael Fall 

Space Telescope Science Institute, 3100 San Martin Drive, Baltimore, 
^ ; MB 21218, USA 

o \ 

Abstract. This article outlines a method to interpret the extragalactic 
background light in terms of the large-scale, average properties of galax- 
S . ies, including the comoving densities of stars and interstellar gas, metals, 

and dust. These quantities are related by a series of coupled conservation- 
type equations analogous to the equations of galactic chemical evolution. 
This approach enables us to combine observations of the emission and 
absorption in galaxies and thus to relate their average stellar and inter- 
J> ' stellar contents. Applications of the method include predictions of the 

global history of star formation from absorption-line observations and 
corrections to the cosmic UV emissivity for absorption by dust. 
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1. Introduction 



The extragalactic background light (EBL) is an inherently average or "global" 
property of the present-day Universe. It is a direction- averaged, time-integrated 
record of the emission and absorption of photons since recombination (z <^ 1000). 
As such, the EBL is closely related to the global history of star formation and, 
potentially, that of other processes, such as the accretion of matter onto black 
holes. The global history of star formation in turn is closely related to the global 
histories of gas consumption and metal production in galaxies. The purpose of 
this article is to explore some of these relations. In this approach, one largely 
ignores the individuality of galaxies and their internal complexity. The primary 
focus, instead, is on the average properties of galaxies. 

The global approach may be contrasted with others in which the individ- 
uality and diversity of galaxies are regarded as paramount. In the latter, one 
attempts to predict, or at least to describe, the full variety of evolutionary paths 
followed by galaxies of different types, masses, and other properties. This is 
necessary to interpret luminosity functions and number counts. The EBL and 
other global quantities are then obtained by integrating over the these func- 
tions. This approach is more challenging and perhaps more fundamental than 
the global approach in that it ties the interpretation of the EBL to a fairly 
complete understanding of the formation and evolution of galaxies. The global 
approach forfeits some of this insight but is much simpler. That is both its main 
limitation and its main advantage. 
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Figure 1. Evolution of the contents of a large comoving volume of 
the Universe, from the Big Bang to the present, including galaxies 
and the IGM. The wavy lines represent the light emitted by stars, 
AGN, and dust in galaxies; the straight lines represent the light emitted 
by quasars and then partially absorbed or scattered in the ISM of 
foreground galaxies and the IGM. 



2. Global Evolution: Overview 

One of the grand themes to emerge in the last few years is the idea that we may 
be able to determine the global histories of star formation, gas consumption, 
and metal production in galaxies from high redshifts to the present. This idea 
is illustrated in Figure 1, where we sketch the evolution of the contents of a 
large, comoving box. We can conveniently quantify the masses of the different 
constituents of the box by the corresponding mean comoving densities normal- 
ized to the present critical density. We are especially interested in the comoving 
densities of stars, gas (both inside and outside galaxies, i.e., ISM and IGM), 
metals, dust, and black holes: Q s , £li sm , Qigm, ^m, and ^bhi respectively. 
We are also interested in the cosmic emissivity E u , the power radiated per unit 
comoving volume per unit rest-frame frequency v, and the EBL intensity J u , 
the power received per unit solid angle of sky per unit area of detector per unit 
observed frequency v. 

After recombination, our comoving box is filled with neutral, metal-free 
gas with nearly uniform density. Perturbations in this intergalactic medium 
(IGM) eventually condense, probably by gravitational clumping and inflow, into 
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protogalaxies. Stars then form in the resulting interstellar medium (ISM). They 
produce metals and may drive outflows of gas from galaxies. In this way, both 
the ISM and IGM may be enriched with metals. Some of the metals remain in the 
gas phase; others condense into solid dust grains. Black holes form in the nuclei 
of some, perhaps even all, galaxies and, when fueled, can power active galactic 
nuclei (AGN). Some of the radiation emitted by stars and AGN propagates 
freely, while the rest is absorbed and then emitted at longer wavelengths by 
dust. Thus, the radiation we detect from galaxies tells us primarily about their 
star, AGN, and dust contents. The spectra of high-redshift quasars contain 
signatures of the absorption and scattering of radiation by the intervening ISM 
and IGM (absorption lines, reddening, etc). Such observations tell us primarily 
about the composition and comoving densities of the ISM and IGM. 

Exactly how all this happens is not yet known, of course. It should be clear 
from Figure 1 and the commentary above, however, that the constituents of our 
comoving box, including the radiation that propagates through it, are very much 
interrelated. In fact, the corresponding comoving densities must obey a series of 
coupled conservation- type equations. Clearly, £l s , Oj sm , and £li gm must add up 
to fJbaryon; a constant. Similarly, Q, s m , f^ m , and f^ m , the comoving densities 
of metals in stars, the ISM, and the IGM, must add up to VL m . Moreover, 
£l m remains a fixed fraction of tt s on the assumption that the global yield is 
constant and that delayed recycling is negligible (a good approximation in the 
present context). The cosmic emissivity E v depends on the star formation rate 
Q s , black hole fueling rate £lbh, and the amount of reprocessing by dust and 
hence Cld- The EBL intensity is related to the cosmic emissivity by a radiative 
transfer equation; for wavelengths at which the IGM is transparent, this gives 



where to is the present age of the Universe. It is likely that AGN make relatively 
small contributions to E u and J v and can be neglected in a first approximation. 

3. Global Evolution: Equations 

To obtain a set of equations for the evolution of the average stellar and in- 
terstellar contents of galaxies, we begin with the familiar equations of galactic 
chemical evolution. In this approach, the IGM is regarded as a reservoir with 
which a galaxy can exchange material by inflow or outflow. Following standard 
practice, we denote the mass of stars and ISM (gas and dust) in a galaxy by M s 
and M g (= M ism ), the interstellar metallicity by Z (= M™ m /M ism ), the inflow 

or outflow rate by Mj, the metallicity of inflowing or outflowing material by 
Zf, and the nucleosynthetic yield of heavy elements, averaged over the stellar 
initial mass function, by y. In the approximation of instantaneous recycling and 
Z <C 1, these quantities are related by 




■to 



(1) 



-(M g + M s ) = M f , 



(2) 



d_ 

It 



(ZM g ) = y-M a -Z-M 8 + Z f M f . 



(3) 
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Equation (3) gives the total rate of change of the mass of metals in the ISM; 
on the right-hand side, the first term is the gain by stellar production and 
ejection (proportional to the rate of star formation since most of the metals are 
synthesized in short-lived stars), the second term is the loss by removal of the 
ISM (including its metals) by forming stars, while the third term is the gain or 
loss by exchange with the IGM (see Tinsley 1980 or Pagel 1997 for complete 
explanations and derivations). 

If we now sum equations (2) and (3) over all galaxies in a large comoving 
volume and then divide by that volume and the present critical density, we 
obtain what are sometimes called the equations of cosmic chemical evolution: 

±(n g + n a ) = n f , (4) 
j t (zn g ) = y e j t Q s -Zj t n s + z f n f . (5) 

Here, the effective yield y e depends on the nucleosynthetic yield y and the spread 
in the interstellar metallicities of galaxies weighted by their star formation rates. 
In the hypothetical case that all galaxies have the same interstellar metallicity, 
y e is the same as y. In general, the two yields will differ, but this is not a 
serious problem because one can adjust y e so that the mean metallicity in the 
models agrees, at some particular time, such as the present, with the observed 
value. This is the same approach taken in nearly all models of galactic chemical 
evolution. Indeed, if one were modeling the evolution of a chemically inhomo- 
geneous galaxy, the same issues would arise in the derivation and solutions of 
equations (2) and (3); the parameter y would then have to be regarded as an 
effective yield, not equal to the nucleosynthetic yield. 



4. Absorption-Line Systems 

The average interstellar properties of galaxies can be determined from the statis- 
tics of quasar absorption-line systems as follows. Let f(N x ,z) be the column 
density distribution of particles of any type x that absorb or scatter light. These 
might, for example, be hydrogen atoms (x = HI), metal ions (x = m), or dust 
grains (x = d). By definition, Hq(1 + z) 3 \dt/dz\f(N x , z)dN x dz is the mean 
number of absorption-line systems with column densities of x between N x and 
N x + dN x and redshifts between z and z + dz along the lines of sight to randomly 
selected background quasars. These lines of sight are very narrow (the projected 
size of the continuum emitting regions of quasars, less than a light year across) 
and pierce the absorption-line systems at random angles and impact parameters. 
One can show that the mean comoving density of x is given by 

= ^TiT^ / N x f(N x ,z)dN x , (6) 
6cHq Jo 

where m x is the mass of a single particle (atom, ion, or grain). Equation (6) 
plays a central role in this subject. It enables us to estimate the mean comov- 
ing densities of many quantities of interest without knowing anything about the 
structure of the absorption-line systems. In particular, we do not need to know 
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their sizes or shapes, whether they are smooth or clumpy, and so forth. A corol- 
lary of equation (6) is that the global interstellar metallicity, Z = ^l 1 ^ 71 /fli 8m , 
is given simply by an average over the metallicities of individual absorption-line 
systems weighted by their total column densities. 

The absorption-line systems of most interest in the present context are 
the damped Lya (DLA) systems. It is widely believed that they constitute 
the ISM of galaxies and protogalaxies and hence are the principal sites of star 
formation in the Universe (Wolfe et al. 1986). There are excellent reasons to 
adopt this as a working hypothesis. First, the DLA systems have, by definition, 
A^hi ^ 10 20 cm" 2 , and this, at least at low redshifts, is near or slightly below the 
threshold for star formation (Kennicutt 1989). Second, the DLA systems contain 
at least 80% of the HI in the Universe and appear to be mostly neutral. The 
other absorption-line systems, those with Ahi <^ 10 20 cm" 2 , probably contain 
more gas in total than the DLA systems, but this must be diffuse and mostly 
ionized. In the following, we regard non-DLA systems as belonging to the IGM, 
even though some of them might actually be associated with the outer, tenuous 
parts of galaxies. This distinction — between the mostly-neutral ISM, where stars 
form, and the mostly-ionized IGM, where they do not — is certainly valid at the 
present epoch (Zwaan et al. 1997). Thus, the DLA systems are often referred to 
as DLA galaxies. The precise nature of these objects — whether they are large or 
small, disk or spheroid — remains to be determined, probably by direct imaging 
(Le Brun et al. 1997). However, as we have already emphasized, the global 
properties derived from equation (6) are not affected by these ambiguities. 

The sample of known DLA galaxies now includes about 100 objects. Most 
of these have redshifts in the range 1.6 <; z <J 4; only a few are known at 
z ^ 1.6 (which require observations from space) and at z ;> 4 (which require 
bright quasars at higher redshifts). We now have measurements of the column 
densities of neutral atomic hydrogen (HI) in all these systems and of the metals 
and dust in many of them (for HI, see Lanzetta, Wolfe, & Turnshek 1995; Wolfe 
et al. 1995; Storrie-Lombardi, McMahon, & Irwin 1996; Rao & Turnshek 2000; 
Storrie-Lombardi & Wolfe 2000; for metals, see Pettini et al. 1994, 1997a, 1999; 
Lu et al. 1996; Boisse et al. 1998; Prochaska & Wolfe 2000; for dust, see Pei, 
Fall, & Bechtold 1991; Kulkarni, Fall, & Truran 1997; Pettini et al. 1997b; 
Welty et al. 1997; Vladilo 1998). It is believed that the DLA galaxies have low 
fractions of ionized hydrogen (HII) because the HI layers are opaque to ionizing 
radiation. In the few DLA galaxies that have been searched, the abundance 
of molecular hydrogen (H2) lies below the value in the Milky Way (Levshakov 
et al. 1992; Ge & Bechtold 1997). One factor that potentially biases all these 
observations is the presence of dust in DLA galaxies. Quasars behind dusty DLA 
galaxies will be obscured, thus reducing the chances that they will be included 
in optically selected samples (Ostriker Sz Heisler 1984; Fall & Pei 1993; Boisse 
et al. 1998). Gravitational lensing has the opposite effect, but this appears 
to be negligible in the existing samples (Le Brun et al. 1997; Perna, Loeb, & 
Bartelmann 1997; Smette, Claeskens, & Surdej 1997). Finally, we emphasize 
that estimates of the comoving densities of HI, metals, and dust are dominated 
by relatively few systems — those with the highest column densities. As a result, 
they are less certain than is often recognized. 
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5. PF Models 

Global evolution models of the type described above have been constructed by 
Lanzetta et al. (1995), Pei & Fall (1995, hereafter PF), Malaney & Chaboyer 
(1996), and Pei, Fall, & Hauser (1999, hereafter PFH). Some related issues are 
discussed by Songaila, Cowie, & Lilly (1990), Timmes, Lauroesch, & Truran 
(1995), Fall, Chariot, & Pei (1996), Calzetti & Heckman (1999), and Harwit 
(1999). In the following, we present highlights from the PF and PFH models. 
At the time the PF models were published, most of the successful surveys for 
emission from galaxies were confined to z <J 0.3, and nothing was known about 

the global history of star formation, tl s {z). Indeed, one of the goals of the PF 
models was to predict £l s (z) from absorption- line observations of DLA galaxies, 
particularly Qm(z) and Z{z). In 1994-95, the comoving density of HI was 
found to be a smoothly increasing function of redshift from z = to z 3, 
while the mean metallicity was known at only two redshifts: Z ~ Z & at z = 
and Z « 0.1Z & at z « 2. These observations — the first pertaining to gas 
consumption, the second to metal production — strongly suggested that most 
stars in the Universe formed after z « 2. The motivation behind the global 
evolution models was to provide a framework to interpret such observations 
more quantitatively. 

The PF models are based on the following simplifying assumptions: 1. HI is 
a tracer of the total interstellar content of galaxies, both atomic and molecular 
(£l g oc f^Hi)- 2. The dust content of galaxies is proportional to the metal content, 
both locally and globally (A^ oc ZN g and oc ZQ g ). 3. The interaction 
between galaxies and the IGM takes one of three idealized forms: closed box 
(Qf = 0), metal-free inflow (Clf = +v£l s , Zf = 0), or metal-enriched outflow 
{(if = —v(l s , Zf = Z). Assumptions 1 and 2 are consistent with, and indeed 
motivated by, the observations summarized in the previous section. Assumption 
2 couples the internal absorption and the obscuration of background quasars 
to the chemical evolution of galaxies; this enables self-consistent corrections for 
biases in the the census of DLA galaxies and hence in the estimates of £lm{z)- 
Assumption 3 is the cosmological analog of the standard treatment of inflow 
and outflow in models of galactic chemical evolution, beginning with the work 
of Larson (1972) and Hartwick (1976). While one might question the precise 
validity of these assumptions, they should be realistic enough to reveal a first 
glimpse of the global evolution of galaxies. 

The PF models were designed to reproduce (as input) the observed evolution 
in Qni{z). As a bonus, they also matched (as output) the available estimates 
of Z{z) without any fine tuning of parameters. Figure 2 shows the predicted 
evolution of the comoving rate of metal production pz{z) i n the standard PF 
models. (The results displayed in this and subsequent figures are for Hq = 
50 kms _1 Mpc _1 , qo = 0.5, and A = 0.) This is just another way of expressing 
the star formation rate: pz = u{3Hq /8ttG)(I s . Figure 2 also shows subsequent 
emission-based estimates of pz{z) from several surveys, including the Canada- 
France Redshift Survey and the Hubble Deep Field (Gallego et al. 1995; Lilly 
et al. 1996; Madau et al. 1996, 1998; Connolly et al. 1997). These were derived 
from rest-frame Ha and UV emissivities and the approximate proportionality 
between UV emission and metal production in stellar populations. Evidently, the 
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Figure 2. Comoving rate of metal production as a function of red- 
shift. The curves are from the PF global evolution models with the 
standard parameters (upper curve with inflow, lower curve with closed- 
box or outflow). The data points correspond to cosmic Ha and UV 
emissivities from Gallego et al. (1995), Lilly et al. (1996), Madau et 
al. (1996, 1998), and Connolly et al. (1997). 

predicted and observed rates are in broad agreement (within factors of about 
two). This is remarkable because the PF models were constructed only with 
absorption-line observations in mind, before the emission-based estimates of 
pz became available. The PF models, when combined with stellar population 
synthesis models, also predicted a far-IR/sub-mm background (Fall et al. 1996) 
that was consistent with early empirical limits and estimates from the DIRBE 
and FIRAS experiments on COBE (Hauser 1996; Puget et al. 1996). 

6. PFH Models 

The PFH models were designed to improve upon the PF models by addressing 
several subsequent developments in this rapidly advancing field. First, it was 
realized that estimates of the comoving star formation and metal production 
rates from rest-frame UV emissivities were potentially biased by absorption by 
dust. The proposed corrections to the estimates of pz in Figure 2 ranged from 
nearly nothing to more than an order of magnitude. The large uncertainty in 
these corrections stemmed from the difficulty of determining the dust content 
and wavelength dependence of the effective absorption in high-redshift galaxies. 
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Figure 3. Extragalactic background intensity J u times frequency u 
as a function of wavelength A. The curves are from the PFH global 
evolution models. References to the observations are given in the PFH 
paper. 



Second, it was found from a large-area, ground-based survey that even the uncor- 
rected UV emissivity at z ;> 3 was probably about twice as high as had initially 
been estimated from the small-area Hubble Deep Field (Steidel et al. 1999). 
Third, the limits on, and measurements of, the EBL intensity have improved, 
especially in the far-IR/sub-mm spectral region (Fixsen et al. 1998; Hauser et al. 
1998; Schelgel, Finkbeiner, & Davis 1998). These are shown in Figure 3. It now 
appears that there is at least as much energy in the long-wavelength hump of 
the EBL spectrum (A ;> 10 ^m) as in the short-wavelength hump (A <^ 10 /xm). 
This, of course, provides a valuable constraint on the amount of UV/optical 
emission that has been absorbed and reradiated by dust. 

In the PFH models, the observed emissivity at rest-frame UV wavelengths 
E u {z) is treated as an input function, along with Clm(z). The emissivity at 
other wavelengths is computed from stellar population synthesis models and a 
template spectrum of dust emission derived from the local emissivity at several 
IR wavelengths. The absorption per mass of dust, assumed proportional to the 
mass of interstellar metals, is then adjusted to satisfy the empirical constraints 
on the EBL intensity. This procedure closes the equations of cosmic chemical 
evolution without invoking a specific relation between the star formation rate 
and the inflow or outflow rate, as in the PF models. Indeed, in the PFH models, 
both Cl s (z) and Clf(z) are independent output solutions of the equations. The 
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Figure 4. Comoving rate of star formation as a function of redshift. 
The open circles represent the observed rate, based on rest-frame UV 
emissivities, from Steidel et al. (1999), while the filled circles represent 
the rate corrected for absorption by dust in the PFH models. 



net effect of this procedure is to determine the correction for absorption by dust 
between the observed and true UV emissivities, i.e., between the observed and 
true star formation rates, as a function of redshift, in a way that is consistent 
with the production of metals and dust and the consumption or replenishment of 
gas by a combination of star formation and inflow or outflow. The results are also 
guaranteed, by construction, to be consistent with a wide variety of observations, 
including the UV emissivity and comoving density of HI as functions of redshift 
and the relative amounts of energy in the long- and short-wavelength humps of 
the EBL spectrum. 

As Figure 3 shows, the EBL intensity in the models does indeed match the 
observations. (In Figures 3-5, the short-dash, solid, and long-dash lines indicate 
the range of solutions permitted by uncertainties in J„ at sub-mm wavelengths.) 
Figure 4 shows the evolution of the comoving rate of star formation Cl s (z) before 
and after the corrections for absorption by dust. Evidently, these corrections are 
factors of 2-4 at most redshifts. The star formation rates in the PFH models have 
the same qualitative behavior as those in the PF models, including the rapid rise 
from z = to z = 1-2, but the decline from z = 1-2 to z 4 is now shallower 
or possibly non-existent. Figure 5 shows the predicted evolution of the mean 
interstellar metallicity and the comoving density of interstellar metals, Z(z) and 
fi^ m (z). The former increases monotonically with decreasing redshift, while the 
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Figure 5. Mean interstellar metallicity in units of the solar value 
(upper panel) and comoving density of interstellar metals (lower panel) 
as functions of redshift. The curves are from the PFH global evolution 
models. The data point at z = 0.8 is from Boisse et al. (1998), while 
those at z > 1.7 are from Pettini et al. (1997a). 

latter first increases and then decreases as the interstellar medium is consumed 
by star formation. The effective yield in the models has been adjusted to give 
Z = Zq at z = 0, the approximate mean interstellar metallicity in present-day 
galaxies (averaged over internal gradients and the luminosity function of galax- 
ies). At higher redshifts, the mean metallicities in the models are marginally 
consistent with those in the DLA galaxies, although the observed enrichment 
appears to be slower than predicted. This may in part be a consequence of the 
bias against selecting dusty and hence metal-rich DLA galaxies, an effect that 
tends to make the observed Z(z) relation artificially flat. 

7. Conclusions 

The main conclusion to be drawn from this article is that we now have the 
apparatus needed to interpret simultaneously a wide variety of observations 
pertaining to the average properties of galaxies, including the comoving densities 
of stars and interstellar gas, metals, and dust. An attraction of the global 
evolution models is that they relate observations that would otherwise appear 
disparate, namely, those of the light galaxies emit and those of the light they 
absorb (from background quasars). Moreover, this approach is perhaps the 
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simplest and most natural one in which to interpret the EBL. Simplicity is 
achieved by focusing on the average properties of galaxies and, for the most 
part, ignoring their individuality and internal complexity. 

The global evolution models have now been developed to the point that 
the accuracy of the results is limited more by observational than by theoretical 
uncertainties. The PFH models incorporated the best data available in 1998. 
Since then, there have been debates about the evolution of the comoving density 
of HI and the mean metallicity in DLA galaxies. The evolution of the UV 
emissivity, and hence the comoving rate of star formation, is probably equally 
uncertain, even before the corrections for absorption by dust. The problem is 
that all these observational relations are based on relatively small samples and 
are affected by a variety of selection biases (e.g., dust for absorption, surface 
brightness for emission). Thus, it should be a high priority in the next few years 
to improve the observations, by increasing the sizes of the samples (e.g., by 
factors of 10 or more) and by reducing the selection biases (e.g., by using radio 
rather than optically selected quasars). It should then be possible to repeat 
much of the analysis described here with greater confidence in the accuracy of 
the results. 
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